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Abstract 

We find a first-order partial differential equation whose solutions are all ultralo- 
cal scalar combinations of gravitational constraints with Abelian Poisson brackets 
between themselves. This is a generalisation of the Kuchaf idea of finding alterna- 
tive constraints for canonical gravity. The new scalars may be used in place of the 
hamiltonian constraint of general relativity and, together with the usual momentum 
constraints, replace the Dirac algebra for pure gravity with a true Lie algebra: the 
semidirect product of the Abelian algebra of the new constraint combinations with 
the algebra of spatial diffeomorphisms. 
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1 Introduction 



In the hamiltonian formulation of general relativity one assumes that the spacetime manifold M , 
with metric 7q,^(X), can be foliated as S x i? where E is a three-dimensional manifold and the 
real line provides a global time coordinate.^ The physical information of the four-dimensional 
theory is then contained in the intrinsic metric gij{x) on S — the pullback of 7q/3(X) on E — and 
its conjugate momentum p^^{x). 

On the foliated manifold, the Einstein equations of general relativity become ten equations 
that the canonical data igij{x),p'''^ (x)) must satisfy. Six of these are first-order dynamical equa- 
tions and describe how the spatial metric and its conjugate momentum change with time; the 
remaining four are constraints that {gij{x),p^^{x)) must obey at any point x on S and reflect the 
symmetries of the four-dimensional theory. The momentum constraints 7ii{x) generate spatial 
diffeomorphisms on a slice S, and the hamiltonian constraint 'H_l(x) propagates E in the direc- 
tion normal to the hypersurface. As functions of gij{x), p^^{x) and the determinant of the spatial 
metric, g{x)^ the gravitational constraints are|jl], |2|: 

n^ix-.g^p] = G,jki{x;gV\x)p''\x)-g"^{x)''R{x-g] (1) 

where Gijki{x;g\ = -g'^{x){gikgji + gugjk - gijgki) 
ni{x;g,p] = -2Djp^^{x). (2) 

The constraint system, 

n±{x) = o = ni{x), (3) 

satisfies the Dirac algebra 

{ndx),n^{x')} = g'^{x)ni{x)6j{x,x')-{x^x') (4) 

{n±{x),ni{x')} = n±^,{x)6{x,x')+n±{x)6,^{x,x') (5) 
{ni{x),nj{x')} = nj{x)5,^{x,x')-{ix^ jx'). (6) 

The Dirac algebra is often thought of as the algebra of spacetime diffeomorphisms, LDiffM, 
"projected" on the foliation T, x R. Equation (^) is the algebra of spatial diffeomorphisms, 
generated by Tli{x). The Poisson bracket (|5|) shows that 7{±{x) transforms as a scalar density 
of weight 1 under spatial diffeomorphisms. The first equation, (Q), however, means that two 
infinitesimal normal deformations on E, performed in arbitrary order, end on the same final 
hypersurface but not on the same point on that hyper surface]^. The fact that the right hand 
side of this Poisson bracket involves the metric explicitly is a source of problems in any attempt 
to use the Dirac algebra in a canonical quantisation of gravity. The Dirac algebra is closed, in 
the sense that the Poisson brackets between constraints are constraints themselves, however, the 
structure coefficients in equation (^) depend on the canonical variables. In a quantum version of 

^ The notation used in this paper is: events on M are denoted by X", while spatial points on S are 
labelled by x^. Indices a,/5, ... are spacetime indices running from to 3; ... are indices on E that run 
from 1 to 3. A standard notation is used to distinguish functions from functionals, namely, round brackets, 
( ), indicate that the object is a function of the arguments in the brackets, square brackets, [ ], denote a 
functional and an object with a mixed bracket, ( ; ], is a function of the arguments that appear on the left 
and a functional of those on the right. Finally, the covariant derivative on (S, (?) is written as Di. 
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the theory where, presumably, quantum counterparts of the classical gravitational data and the 
constraints will appear, the Poisson bracket algebra will go over to a commutator algebra which 
is not a Lie algebra. 

Recently, Brown and Kuchafj^ obtained some surprising results which provide a promising 
proposal on the problem of the algebra of canonical gravity. They studied a spacetime filled with 
incoherent dust; in such a system the coupling of the metric to matter introduces into spacetime 
a privileged dynamical frame and time foliation. The system has a very interesting feature: the 
scalar constraint for the combined system, which in was denoted by Hi^{x), can be split into 
two separate parts, one for matter and one for gravity, and cast in the form 

H^{x) := P{x) + h{x; g^„p'^]=0, (7) 

where P{x) is the momentum conjugate to the "dust time" variable T(x), and h(x) is the following 
scalar combination of the gravitational constraints: 

h{xf = G{x) := nl{x) - g'^rLi{x)rLj{x). (8) 

The truly remarkable result of ^ was that G{x) and h{x), scalar functions of gravitational 
variables only, have strongly vanishing Poisson brackets among themselves. One then has a 
possible alternative set of constraints for pure gravity, 

G(x) =0 = (9) 

which generate the true Lie algebra, 

{G{x),G{x')} = 0, (10) 
{G{x),ni{x')} = G,i{x)5{x,x') + 2G{x)5^i{x,x'), (11) 

{ni{x),nj{x')] = njix)64x,x')-iix^jx'). (12) 



This algebra is the semidirect product of the Abelian algebra generated by G{x) (equation [LOD , 
and the algebra of spatial diffeomorphisms LDiff S generated by T~Li{x) (equation . The Poisson 
bracket (|llj) is the transformation of G{x) as a weight 2 scalar density under Diff S. 

Previous investigations on the use of matter time and reference fiuids in canonical gravity]^, |^ 
led to the "phenomenological" approach to the problem of time in canonical gravity of ||5|. The 
significance of the new constraints (^), however, and the role they could play in pure gravity was 
left to be investigated at some later stage. That dust was not, in itself, essential in the process 
was made clear when Kuchaf and Romano[^ found a simpler way of producing results similar 
to those of ref. ||5|. They coupled gravity to a massless scalar field, and the scalar density they 
obtained is the following combination of gravitational constraints 

K±{x)=g"2{x)(^-n^{x)±^G{x)^. (13) 

As in the case of G(x), the weight -two scalar densities A-i-(x) have strongly vanishing Poisson 
brackets with themselves. They can be used to rewrite the constraint system for pure gravity as 

A±(x) = = W,(x) (14) 
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and, again, if A-i-(x) can be taken to replace a true Lie algebra for vacuum gravity has 

been found. 

Although coupling to matter was an essential step for ||5| and to find the explicit form 
of the new scalar constraint, the resulting combinations, G{x) and A(x), involve gravitational 
variables only and do not depend on the matter variables to which the gravitational field was 
originally coupled. The question that naturally arises therefore is whether a transformation of 
the scalar constraint to a form that generates a true Lie algebra can be made within the context 
of pure gravity only. 

In this paper we take the first step towards answering this question by finding the full set 
of scalar densities, combinations of the original pure gravity constraints, that have the crucial 
property of generating a Lie algebra (with vanishing Poisson brackets between themselves). We 
find that the calculation can be carried out by working directly with the Dirac algebra, without 
any reference to matter couplings. Thus, we are not restricted by any need to include matter 
fields and there is greater freedom in using the new algebra in pure gravity. A complementary 
view of these combinations of pure gravity constraints is provided in a paper by Kouletsis|^], to 
appear shortly. In this work the constraint combinations are shown to arise naturally in a system 
of pure gravity with non-derivative coupling to an action of a single scalar field and two arbitrary 
functions of a Lagrange multiplier. 

2 Constraint combinations are generated by a 
differential equation. 

From the work of Q and ^ , it seems clear that there might be more combinations of gravitational 
constraints obeying the algebra (p!o|-|T^). We want to find them all, and furthermore, as explained 
in the introduction, we wish to investigate the new algebra without any reference to a matter 
coupling. We therefore choose to work directly with the Dirac algebra of pure gravity. 

A particularly interesting feature of the constraints G{x) and A±(x) is that they are scalar 
densities of weight 2, unlike the usual hamiltonian constraint 7i^{x) which is a weight 1 scalar 
density. An important question therefore is whether a weight of 2 is in some way a natural choice 
for a scalar constraint for pure gravity that has Abelian Poisson brackets with itself. 

For generality then, we assume that the candidate new constraint we are looking for is an 
ultralocal scalar density /C(x) of arbitrary weight oj. Such a }C{x) can only be a function of the 
scalars 7i±, {g^^TiiTCj), and the determinant of the spatial metric g (there are no derivatives of 
the gravitational constraints). Together with the usual momentum constraint, /C(x) is assumed 
to satisfy the Lie algebra 

{/C(x),/C(x')} = 0, (15) 
{IC{x),7ii{x')} = IC^i{x)5{x,x') + LjlC{x)5^i{x,x'), (16) 

{niix),njix')} = nj{x)6,i{x,x')-{ix^jx'). (i7) 

The bracket (|l^), consistent with our original assumptions, means simply that }C{x) transforms 
as a scalar of weight oj under DiffS. The dependence of /C(x) on 7i±{x) and TCi{x) lies in the 
first bracket (ITsI), which we shall now evaluate. 

The calculations that need to be done are considerably simplified if we define h and / as the 
two obvious scalar densities of weight zero that can be formed from the gravitational constraints 
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T-L±_ and Hi: 

h:=g-hni_ and f := {g~^)g'minj. (18) 

Using h and /, one can define another weight zero scalar density K{x), related to the weight uj 
density /C(x) by 

}C{x;h,f,g]=g^/\x)K{x;h,f]. (19) 

It is easy to construct a IC{x) by premultiplying K{x) by the g^^'^{x) factor, so for simplicity, our 
discussion is centred around K[x). One only needs to remember that K{x) itself does not satisfy 
the Abelian Poisson bracket (p!^). Its dependence on /i, /, however, is the same as that of /C(x). 

Imposing the central requirement that /C satisfies the Poisson bracket (|l5|), we can now find, 
by simply working out the Poisson brackets, which (weight zero densities) K correspond to the 
(weight u densities) /C that generate a true Lie algebra . Our calculation is a generalised analogue 
of the calculation of Kuchar^ and Romano, so we concentrate on the steps needed to adapt their 
calculation for arbitrary- weight densities. 



By the Leibniz rule, equation (15) gives 

{/C,/C'} = g^'^g^'\K,K'} + [g'^'^K{9^'\K'} - (x ^ x')) , (20) 

where, to avoid cluttering of symbols, primes have been used on quantities that have x' as their 
spatial argument. The Poisson brackets above can be expanded to 

{9-/',K'} = ^{9^^^h'} + ^{g-/M'}, (22) 

and therefore, one needs to calculate Poisson brackets involving h, f and g^^'^. To handle these 
brackets, the identity 

{g''/\f'} = cog'^{gKf'}, (23) 

is used, while the bracket between g^^'^ and h is proportional to 5{x, x') and therefore its con- 
tribution cancels with the {x x') terms that appear in equation (^). For the brackets of the 
weight zero gravitational variables h and /, it is straightforward, using the Jacobi identity and 
the standard Dirac relations (§)-(§), to verify that they satisfy the algebra 

{h,h'] = g-h''Hn±,n'^}, (24) 
{/,/'} = -ig-h'-^^f{n^,n'^}, (25) 

{hJ'} = (x6{x,x'). (26) 

where oc 6{x,x') stands for terms proportional to 5{x,x') which, as always, cancel with the 
(x ^ x') terms. One now has all the necessary information to extend the procedure of Kuchaf 
and Romano[P for /C(x). 
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The requirements on the scalar density IC{x) are minimal. Thus one would expect that 
the rather lengthy calculation of the Poisson brackets between the /C's would lead to a fairly 
cumbersome result. It is surprising to find it is instead quite simple: 

{n±,n'^}. (27) 

All functions K that turn the quantity inside the square brackets in ( p7| ) into zero correspond, via 
equation (p!9|), to candidate new constraints fC. The family of these scalar functions are therefore 
the solutions of the nonlinear, first-order partial differential equation for K 




This unexpectedly compact equation is the key result because any of its solutions can be 
directly used to give a Lie algebra for pure gravity. Also, if by some other method one has 
already obtained a combination of gravitational constraints which is suspected to satisfy a Lie 
algebra, it is enough (and far easier) to instead check that it satisfies the differential equation. 
For example, it can be easily verified that the differential equation ( p8|) is satisfied by the scalars 
that were discovered in references |5| and (in their respective weight zero form), 

K = g-^G = h^ - f and K = g~^A± = -h± ^ h? - f , (29) 
for weight lo = 2. It is also satisfied by the weight uj = \ scalar yh^ — f that appeared in ref. 



{/C,/C'} = <7' 



3 The solution of the differential equation 

The general solution to the first-order nonlinear partial differential equation (^) of the two 
independent variables h and /, for 7^ 0, is expected to be an expression involving Ti.^, 7ii and 
an arbitrary function of one parameter. 

A change of variables can turn (^) into a more manageable equation. If we set 

C = lnK (30) 

(we assume that C in the above definition is well-defined) and then divide the differential equation 
( p^ through by K^, we can rearrange and rewrite it as a partial differential equation for C: 

The advantage of this rearrangement is that this equation does not involve C explicitly. Be- 
cause the variables h, f are independent, the differential of the function C{h,f) is given by the 
expression 



dC(hJ) 
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In equation there are now no terms involving both variables h and /, and thus it can be 



"split" into two parts, for dC/df and dC/dh; namely, it holds |11| that any solution to the 
coupled set 



\df) ~ 2df 




= (33) 

for some constant a, will also satisfy (|3lD. These two, single- variable, equations can be readily 
solved as 

— = - (^a; ± v/JTTe^ j and ^ = ±2a (34) 

(note that the it signs appearing in the above pair of equations are independent) and therefore 
the differential of C{h,f ) in equation (5^) can be integrated to 



C{h,f,a,l3) = jlnf±^ 



2a/lj2 + 16a2/ + a;ln 



yjufi + \^C? f — to 

^uji'^ + 16q2/ + W 



± 2a/i + /? (35) 



(where, again, the two possible choices of it sign are independent). This is a complete integral of 
the differential equation for C{h,f); namely, a solution of the equation involving two arbitrary 
independent constants a, [3 which come from the two integrations in equation (|32[). The complete 
integral then is a two-parameter family of surfaces, since a different surface is obtained for each 
choice of a and /3. 

The general solution should also describe this family of surfaces. The difference is that, 
instead of depending on two arbitrary parameters, the general solution should contain an arbitrary 
function of one parameter. That is, we suppose that 

/3 = (/>(«) (36) 

so that C{h, /, a, (^{a)) describes a one-parameter family of solutions. One can think of (/'(a) as the 
two- variable analogue of the integration constant of the familiar one- variable partial differential 
equation. 

One can go further by bringing in a basic theorem in the theory of first-order partial dif- 



ferential equations |11], that the envelope of any family of solutions of a first order equation, 
depending on some parameter, is again a solution. The envelope of the one-parameter family 
of solutions C(/i, /, a, 0(a)) is essentially a surface tangent to the family of surfaces that are the 
solutions. It should then be given by C{h, f,a,(j){a)), together with its differential with respect 
to a, namely, 

— = ± ±2h + cp{a) = 0, (37) 

where a prime denotes differentiation with respect to a. One now only has to eliminate a from 
the above equation and C{h, f,a,(j)(a)) to have a single expression for the envelope. This will 
also be the general solution involving the arbitrary function (/)(a)[^]. If we return to the original 
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equation for K by exponentiating C in equation (35), we can write out the general solution to 
the differential equation (]28| ) directly. It is given by: 



i^±(/i,/,0(a)) 



h±\<l,\a)\±J{h±]^<l,'{a) 



2 



X exp ( 0(a) T ^^^i^==\ (38) 




(where the it signs in front of </>(«) are independent of those in front of the square roots) together 
with equation (|37[). Note that, by a trivial rewriting of the exponential factor, the whole solution 
can be raised to the power of (aj/2), and hence the constraint combination /C(x), is simplest but 
non-trivial when uo = 2. This is consistent with the fact that the "physical" matter couplings 
used by Kuchaf et al, also produced weight 2 densities. 

Therefore, for each choice of the function (p{a), one solves (^) algebraically to obtain 0(a) as 
a function of h and /, and then substitutes the result in the general solution (pSj). For example, 
the Brown-Kuchaf solution is reproduced when (f){a) = — In a — 1 and uj = 2 (and the minus signs 
in front of the square roots), so that (j)'{a) = — ^- Then, solving ( |37D gives 

(^ = -j;TZrj^ (39) 
which implies (j)'{a) = h — t, and which, when inserted in (^), produces K = ■ii{h? — /). 



Obtaining the Kuchaf-Romano solution is also straightforward, it only requires setting 0(a) = 0. 
With the same method one can also generate new solutions, in fact an infinite number of them, 
for all possible choices of 0(a). New solutions, and also their relationship to scalar field actions, 
have been investigated by Kouletsis in [^.^ Generally, the only simple ones are those that have 
been already found by Kuchaf et al. 

The weight u; = case has to be treated separately since then equation ([2^) becomes a 
homogeneous partial differential equation. 




1 (dK'^V 



dh 



0, (42) 



whose solution is of the form = oc (/i it ^/f) + constant. 

Summarising, we have shown that for each choice of the function 0(a) in the pair of equations 
(|37| ) and (|38|), a combination of the gravitational constraints h and / is generated. The solution, 
via equation (|l^), provides the Abelian weight oj scalar density /C(x) which can be used in the 
construction of a true Lie algebra for pure gravity. 

^ For example, the function 0(a) = Alna, where A is an arbitrary constant, leads to the new 
solution lH^: 

where 



//12 I 1 — ^ 2 

Ri = ^h^ + (A2 - 1)/ and R^^l^^h^ + —^f ^-hR,. (41) 
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4 Conclusions 



In this paper we started from observations of the alternative hamiltonian constraints found by 
Kuchaf et al., which generate a true Lie algebra of gravity. We generalised to a much larger family 
of such scalar densities in the hope that knowing their general form would make their origin clearer. 
All previous work on these constraints has been based on introducing some simple form of matter 
into spacetime, either dustp, a scalar field|P or the deWitt action for a perfect fluid[12|. We 
found that it is possible, and in fact simpler, to arrive at these constraint combinations working 
solely with the algebra of constraints. This is also the least restrictive method and enabled us to 
generalise to ultralocal scalar densities of arbitrary weight. Here we saw that, although in principle 
scalar densities of any weight can be employed in the construction of a Lie algebra, the weight 
2 ones are the most simple and natural choice. In the present work, we have not addressed the 
problem of the global equivalence of the scalar densities found to the usual hamiltonian constraint 
of general relativity. We have also not discussed the conditions for the constraint combinations 
obtained to be well-defined. For example, it has already been noticed by Brown and Kuchaf that 
the hamiltonian vector field associated with the new scalars vanishes on the constraint surface. 
The criteria determining the answer to these questions may depend on the context in which the 
new scalars will be applied. 

There are certainly other interesting aspects of these gravitational constraint combinations. 
They are related to matter time in canonical gravity and to coordinate conditions!^. Possibly, 
they can also come out of a canonical transformation of the geometric canonical data that sepa- 
rates the true degrees of freedom (in the spirit of jl^). The possibilities they open in the quantum 
theory are certainly worth investigating and we will discuss them in future work. 
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